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ABSTRACT 

Helix 69 (H69) is a 19-nt stem-loop region from the 
large subunit ribosomal RNA. Three pseudouridine 
(T) modifications clustered in H69 are conserved 
across phylogeny and known to affect ribosome 
function. To explore the effects of 4* on the conform- 
ations of Escherichia coli H69 in solution, nuclear 
magnetic resonance spectroscopy was used to 
reveal the structural differences between H69 
with (¥4*4') and without (UUU) ¥ modifications. 
Comparison of the two structures shows that H69 
TTT has the following unique features: (i) the loop 
region is closed by a Watson-Crick base pair 
between T1911 and A1919, which is potentially 
reinforced by interactions involving T1911N1H and 
(ii) T modifications at loop residues 1915 and 1917 
promote base stacking from 4*1915 to A1918. In 
contrast, the H69 UUU loop region, which lacks T 
modifications, is less organized. Structure modula- 
tion by ¥ leads to alteration in conformational 
behavior of the 5' half of the H69 loop region, 
observed as broadening of C1914 non-exchangeable 
base proton resonances in the H69 TTT nuclear 
magnetic resonance spectra, and plays an important 
biological role in establishing the ribosomal 
intersubunit bridge B2a and mediating translational 
fidelity. 

INTRODUCTION 

RNA molecules can adopt highly folded 3D structures to 
carry out their essential structural and catalytic functions in 
biological systems (1). As enrichment to the four standard 
nucleotides (i.e. A, C, G and U), post-transcriptional modi- 
fications enhance the chemical repertoire of RNA and play 
important roles in 'fine-tuning' local conformations of 
RNA (2,3). Among the > 100 modifications identified 
to date (4), pseudouridine (*f) (Figure la) was the first 
reported and is also the most frequently encountered (5). 



Uridine (Figure lb) is isomerized to * (Figure la) by 
replacing the N-glycosidic bond with a C-glycosidic bond, 
and this covalent structural variation has been shown to 
modulate local conformation and overall activity in tel- 
omerase (6), spliceosomal (7) and transfer (8) RNAs. In 
the Escherichia coli (E. coli) ribosomal RNAs (rRNAs), 
10 ^ modifications are distributed in the functionally im- 
portant regions (e.g. the peptidyl transferase center (PTC) 
and the intersubunit bridge B2a) (9,10), with the latter 
hosting three in a 19-nt-long hairpin segment of the 
23S rRNA named helix 69 (H69). 

Helix 69 shows a high degree of conservation in both 
sequence and secondary structure across phylogeny (11). 
An additional conserved feature of H69 is the existence of 
multiple pseudouridylation sites (E. coli numbering, pos- 
itions 1911, 1915 and 1917; Figure lc and d), which have 
been mapped in E. coli, human and several other organ- 
isms (9,10,12-14). Helix 69 plays a central role, together 
with helix 44 (h44) of the small subunit rRNA (16S and 
16S-like), in establishing the B2a intersubunit bridge, one 
of the most highly conserved regions of the ribosome 
(15-19), and has direct contacts with multiple factors 
during the translation process, as revealed in crystal and 
cryo-EM structures (18,20-26). Together, these facts point 
to the significance of the conserved ^s in folding and 
function of H69 (27,28). 

Mutagenesis studies provided insight into the conserva- 
tion of sequence and/or post-transcriptional modifications 
in H69. Deletion of H69 in bacterial ribosomes leads to the 
disruption of the B2a intersubunit bridge, which in turn 
results in defects in ribosome association and peptide 
release (29). Single-nucleotide mutations in H69 confer 
slow-growth phenotypes with compromised protein synthe- 
sis activity and translational fidelity (30-33). Among the 
H69 loop residues most sensitive to single-nucleotide mu- 
tations, *1915, A1916 and *1917 show reduced growth 
rates and interrupted ribosome assembly (33). It was 
revealed in the post-transcriptional modification mapping 
of H69 that the mutation A1916G reduces isomerization of 
U1917 to *1917, and A1916U abolishes the modification 
activity completely (34). Therefore, the slow-growth pheno- 
type observed with A1916G is partially attributed to 
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Figure 1. Sequence and modification sites of H69 from E. coli 23S 
rRNA are shown, (a) A pseudouridine (*) has a C5-C1' glycosidic 
bond, (b) A uridine residue contains an Nl-CT glycosidic bond. The 
numbering for the imino protons of U and * is given, (c) The RNA 
construct without pseudouridylation, UUU, is given with the E. coli 
numbering, (d) The *** construct shows pseudouridine (*) residues at 
positions 1911, 1915 and 1917. 



reduced ribosome assembly due to lowered 4»1917 modifi- 
cation levels. A growth advantage with pseudouridylation 
at positions 2258 and 2260 in Saccharomyces cerevisiae 
(corresponding to 1915 and 1917 in E. coli) in the H69 
loop was also reported (35). 

Thermodynamic and chemical-probing studies were 
carried out previously to explore the effects of 
pseudouridylation on H69. It was demonstrated on 
other RNA motifs that * can stabilize RNA folding by 
formation of water-mediated hydrogen-bonding inter- 
actions through ^NIH (36); however, a stabilizing ef- 
fect was only observed at position 1<1911 in H69 model 
oligonucleotide hairpins from E. coli and human. The 
loop-closing *I* contributed —0.6 to — 1 . 1 kcal/mol to 
the AG° 37 of the RNA stem-loop structure (37,38). In 
contrast, 4*1915 and *IT917, individually or collectively, 
showed slight destabilizing effects in the same model 
studies (37,38). Corresponding differences in H69 flexibil- 
ity were observed through SHAPE (39) analysis of 50S 
subunits from E. coli wild-type and RluD - (^-deficient) 
strains; only A1913 and A1918 in the wild-type 23S rRNA 
showed strong reactivity toward the SHAPE reagent, 
whereas all H69 loop residues demonstrated mild reactiv- 
ity in the unmodified RNA (RluD" 23S rRNA) (40). 

To elucidate in more detail the structural effects of * 
modifications on H69 folding and explore correlations 
between the modifications and their biological signifi- 
cance, the solution structures of RNA constructs 
with (*I / <I'*I / ) and without (UUU) pseudouridylations 
(Figure lc and d), representing H69 from E. coli 23S 
rRNA, were examined by using nuclear magnetic reson- 
ance (NMR) spectroscopy. A comparison of the two 
structures reveals that ^s substantially alter the folding 



of the H69 loop region. In UUU, the base moieties of 
all three loop U residues are found to have greater 
solvent accessibility than the corresponding *!> residues in 
which may help with RluD recognition and cataly- 
sis. The 4T91 1 forms a Watson-Crick base pair with 
A1919 and has unique hydrogen-bonding interactions. 
The NMR structure of also shows that ^1915 and 
1*1917 participate in base stacking in the 3' half of the H69 
loop. Together, the three * modifications influence con- 
formational behavior of the 5 1 half of the H69 loop region, 
as shown by line-width broadening of the C1914 base non- 
exchangeable protons, and are suggested to play a role in 
facilitating base flipping of A1913, which is known to 
make important contacts in the B2a intersubunit bridge 
of intact ribosomes (41). 



MATERIALS AND METHODS 

Preparation of H69 RNA oligonucleotides 

Unmodified H69 RNA samples (UUU, 5'-GGCCGUAAC 
UAUAACGGUC-3') were synthesized by in vitro T7 RNA 
polymerase transcription with unlabeled or 13 C, 15 N -labeled 
NTPs, synthetic gel-purified DNA template, and promoter 
sequences (42). Full-length H69 RNA transcripts were 
purified by using denaturing 20% (w/v) preparative poly- 
acrylamide gel electrophoresis and electroelution in 
0.2 x Tris borate + ethylenediaminetetraacetic acid buffer 
with a Schleicher and Schuell® Elutrap. RNAs were 
desalted with Sep-pak® (Waters) reverse-phase chromatog- 
raphy cartridges, and the eluted fractions were pooled and 
lyophilized to a powder. 

Synthetic modified RNA (37) (WW, 5'-GGCCG*AA 
C*A*AACGGUC-3') was purchased from Dharmacon® 
(Thermo Scientific) and subjected to high-performance 
liquid chromatography purification on a Waters Xterra 
MS C18 column. A gradient of acetonitrile from 6.0 to 
7.8% over 24min in 25 mM of triethylammonium acetate, 
pH 6.5, at a flow rate of 3 ml/min was used. The RNA- 
containing fractions were lyophilized and desalted with a 
Sep-Pak column. The molecular masses of the RNA oligo- 
nucleotides were confirmed by using MALDI-TOF mass 
spectrometry. 



Preparation of RNA NMR samples 

Purified H69 oligonucleotides (UUU and ***!>) were 
dissolved in 300 ul of lOmM potassium phosphate, 
50 mM KC1, pH 7.3, and 0.1 mM Na 2 -ethylenediaminete- 
traacetic acid. The samples were lyophilized to a powder. 
Deuterium oxide (D 2 0, 99.96%; Cambridge Isotope Labs) 
was used to exchange the residual H 2 0 twice by lyophiliza- 
tion, and the samples were dissolved to a final volume 
of 300 ul for NMR experiments on non-exchangeable 
protons. A H 2 0/D 2 0 (90%/10%) mixture (300 ul) was 
used for experiments on exchangeable protons. All of the 
NMR samples contained ~1.0mM RNA and a trace 
amount of 3-(trimethylsilyl) propionate (TSP) as an 
internal proton chemical shift reference. 
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NMR spectroscopy 

All NMR experiments were carried out on a Bruker 
Avance 700 MHz spectrometer equipped with an HCN 
cryoprobe and a Varian Mercury 400 MHz equipped 
with a room temperature QXI probe for 3I P NMR experi- 
ments. Spectra for exchangeable proton resonance assign- 
ments and base-pairing identification were acquired at 283 
and 288 K. All other spectra of UUU and samples 
were recorded at 298 and 310 K, respectively. Topspin 2.1 
(Bruker) and Sparky 3.114 (University of California, San 
Francisco, CA, USA) were used for spectral processing 
and analysis. 

The 2D NOESY spectra of the unlabeled H69 RNA 
samples (UUU and W*I>) dissolved in 99.96% D 2 0 
were initially analyzed to assign the resonances of base 
protons (H8/H6) and sugar protons (HI'). Standard 
NMR experiments (e.g. 2D NOESY, 2D DQF-COSY, 
2D 13 C- H HMQC and 3D TOCSY-NOESY) were 
collected on unlabeled samples, whereas 2D 13 C-'H CT- 
HSQC, 3D HCcH-COSY, 3D HCcH-TOCSY, 3D HCCh- 
TOCSY and 3D NOESY-HMQC were run on the 
13 C, 15 N-labeled sample (for UUU only) to assign H/C 
resonances. The ID 31 P and 2D 31 P-'H HETCOR experi- 
ments on unlabeled constructs were used to identify 
chemical shifts of 31 P resonances (UUU and 1 J>V[ / ^I / ). 

Structure calculations 

Interproton distance restraints were derived from 2D 
NOESY spectra as described previously (43). Interproton 
distance restraints of exchangeable protons were only 
applied to residues involved in base pairs in the stem 
region (from G1906-C1924 to G1910-C1920 for UUU 
and from G1906-C1924 to *1911-A1919 for ***) (44). 

The dihedral angles, a, (3, y, §, s, L, and %, of residues in 
the stem region (from G1906-C1924 to G1910-C1920 for 
UUU and from G1906-C1924 to *1911-A1919 for 
were restrained according to the A-form RNA helix 
geometry (-62 ± 10°, 172 ± 10°, 60 ± 10°, 84 ± 10°, 
-160 ± 10°, -71 ± 10°, and -160 ± 10°, respectively). 
Loop residue dihedral angles a and £ were restrained 
(0 ± 120°) to exclude the trans conformation, owing to 
absence of any downfield-shifted P resonances in UUU 
and *W ID 31 P spectra (45,46). Because only weak 
cross-peaks between P-H5' or P-H5" were observed, the 
(3-dihedral angles of loop residues were loosely restrained 
to — 180 ± 90°, which essentially covers the entire steric- 
ally accessible range (47). The y dihedral angles were 
restrained to be 120 ± 120° to exclude gauche' conform- 
ations which are rarely observed. Cross-peaks from 
J-coupling of HI' and H2' in the 2D DQF-COSY were 
used to determine the 8 dihedral angles of loop residues, in 
which an intense cross-peak comparable with the pyrimi- 
dine H6-H5 cross-peaks indicates a CT-endo conform- 
ation, and the 8 dihedral angle was restrained to 
157 ± 40°. If the sugar HF-H2' cross-peak was weak or 
invisible, then the 8 dihedral angle was restrained to 
84 ± 20° corresponding to Ci'-endo conformation; other- 
wise, it was left unrestrained (48,49). The £ dihedral angles 
of loop residues were restrained to be — 120 ± 120° to 
exclude the gauche + conformation. Owing to the absence 



of intense cross-peaks between base protons H8/H6 and 
sugar protons HI' in the 2D NOESY spectra with a short 
mixing time (x m = 60 ms) from UUU and ty^I'ty, the % 
dihedral angles of loop residues were all restrained to 
— 110 ± 110° to exclude a syn conformation. 

The Crystallography and NMR System (CNS) version 
1.2 using a simulated annealing and restrained molecular 
dynamics protocol was used in structure calculations (50). 
Parameterization of * residues was used as previously 
reported (51). An extended structure containing only the 
RNA sequence and covalent linkages of the RNA con- 
struct generated using the CNS 1.2 was subjected to a 
torsion angle molecular dynamics at 20 000 K for 40000 
steps (1 fs/step, 40 ps in total) with molecular dynamics 
force constants of 50 kcal/mol/A 2 and 150kcal/mol/rad 2 
for the distance restraints and dihedral angle restraints, 
respectively. Molecular dynamics pseudopotential force 
constants of 100 kcal/mol/A 2 and 250 kcal/mol/rad 2 were 
then applied to distance and dihedral restraints, respect- 
ively, in the first round (50000 steps in 100 ps) of slow 
cooling down to OK. In the second round (10000 steps in 
35 ps) of slow cooling down from 2000 to OK, a Cartesian 
molecular dynamics simulation was used /using 
pseudopotential force constants of 200 kcal/mol/A 2 for 
the distance restraints and 500 kcal/mol/rad 2 for the 
dihedral restraints). A force constant of 50 kcal/mol/A 2 
was used for the base pair planarity restraints in the 
final energy minimization stage consisting of 800 steps 
for 20 rounds^ (using pseudopotential force constants of 
300 kcal/mol/A 2 for the distance restraints and 700kcal/ 
mol/rad for the dihedral restraints). Qualified output 
structures (i.e. those structures with no distance violation 
>0.3 A or dihedral violation >5°) were chosen for further 
torsion angle molecular dynamics refinement with an 
energy penalty of 5 kcal/mol equivalent to an NOE viola- 
tion of 0.3 A or a dihedral angle violation of 5°. RNA123 
(DNA Software, Inc.) was used for restraint violation 
analysis in NMR structure calculations and refinement 
for all pair-wise root-mean-square deviation (RMSD) cal- 
culations in this work (52). Structure calculation restraints 
and statistics are shown in Table 1 . All the structures were 
visualized with Pymol (The PyMol Molecular Graphics 
System, Schrodinger, LLC). 

RESULTS AND DISCUSSION 

NMR spectroscopy of H69 UUU and T«P«P 

The oligonucleotides used in this study correspond to nu- 
cleotides 1906-1924 in E. coli 23S rRNA, either without 
(UUU) or with (*I / *I'4 / ) pseudouridine modifications at pos- 
itions 1911, 1915 and 1917 (Figure lc and d). In addition 
to a natural abundance UUU sample, a 13 C, 15 N-labeled 
sample of UUU was synthesized by in vitro transcription 
catalyzed by T7 polymerase and used for assignments of 
overlapping resonances. The structure families of the two 
oligonucleotides show comparable quality and conver- 
gence (Table 1). The final solution structures of UUU 
and <I' V I / * are consistent with all of the observed NMR 
data, and support the results of previous structural, bio- 
physical and chemical-probing studies besides providing 
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Table 1. Structure calculation restraints and statistics of UUU 
and *** 



Restraints 


T TT TT T 


I T I 

WWW 


NUh distance restraints 


~)A 1 
-4 1 


__4- 


Intra-residue 


97 


81 


Inter-residue 


144 


143 


Stem constraints 


59 


59 


Base pair planarity 


5 


5 


Dihedral angle restraints 


128 


127 


Statistics 






Violations 






NOE distance violation >0.3A 


0 


0 


Dihedral angle violation >5° 


0 


0 


Average pair-wise RMSD (A) 


1.15 


0.66 


of all atoms (10 structures) 







further information on the structural role played by the *s 
in the formation of the intersubunit bridge B2a. 

Resonances of imino protons from the stem residues of 
both oligonucleotides were observed in the 2D NOESY 
spectra (H 2 0/D 2 0: 90/10%; Figure 2) (53). Cross-peaks of 
the two imino protons from the neighboring base pairs were 
shown to have medium to weak intensity. The chemical 
shifts of G1907N1H and U1923N3H are upheld shifted 
relative to the other imino protons, and an intense cross- 
peak of the two protons exists for both UUU and 
consistent with a wobble G - U mismatch structure. 
Formation of this G*U wobble pair protects the 
G1907N1H and U1923N3H from being exchanged with 
solvent; therefore, the corresponding resonances are 
clearly visible in the 2D NOESY spectra (Figure 2). The 
assignment of the *1911N3H resonance was confirmed 
based on cross-peaks with G1910N1H and A1919H2 
(Supplementary Figure SI a). The strong ^^UNSH to 
A1919H2 cross-peak is consistent with the formation of 
Watson-Crick hydrogen bonding. One upfield-shifted 
imino proton resonance at 10.2 ppm was observed in the 
ID imino proton spectrum of tyl'ty, and this proton is pos- 
itioned in proximity to *1911H6 by an intense cross-peak 
shown in the 2D NOESY (H 2 0/D 2 0: 90/10%) spectrum 
(Supplementary Figure Sib). This resonance is unambigu- 
ously assigned as ^1911N1H and its chemical shift value 
is in good agreement with previous studies (53). No reson- 
ances from other ^ residues in the loop region were 
observed, suggesting that they are exposed to solvent. 

The H8/6 to HI' 'walk' pattern of the UUU stem region in 
the 2D NOESY (D 2 0: 99.96%) spectrum corresponds well 
with that of fl'ty (Supplementary Figure S2a and b). The 
chemical shifts of Hi's of the three * residues in are 
upheld shifted compared with those of the corresponding 
uridine residues in UUU by 0.9, 1.1 and 1.3 ppm, 
respectively, consistent with values reported in other 
^-containing RNAs (6,54-56). The upheld shift of HI' in * 
compared with U can be explained by the fact that the less 
electronegative C5 of ^ replaces Nl of U at the glycosidic 
bond position (57,58). Protons U/*1911H6, A1916H2, 
A1916H1', U/4T917H6 and A1918H2 also show chemical- 
shift differences of >0.1ppm (Supplementary Figure S3) 
between the two constructs. In the 2D NOESY (D 2 0: 
99.96%) spectrum of the line width of the C1914 



H6-H5 cross-peak is significantly broadened compared 
with that of UUU (Supplementary Figure S2c and d). 

Additional dihedral angle restraints were obtained from 
2D DQF-COSY spectra. In the 2D DQF-COSY spectrum 
of UUU, three intense cross-peaks of HF-H2' were 
assigned to C1914, U1915 and U1917 (Supplementary 
Figure S4a), indicating a CT-endo sugar pucker conform- 
ation. In the 2D DQF-COSY of the sample, no Hl' 
H2' cross-peaks are observed. However, the ^I**!** 
spectrum does have one upheld-shifted cross-peak of 
weak intensity corresponding to H6-H1' of 4*1915 
(Supplementary Figure S4b). The appearance of this 
J HH-coupling enables the addition of a dihedral angle 
(H1'-C1'-C5-C6) constraint of 180 ± 45°, as a coplanar 
geometry of H6-C6-C5-C1'-H1' is required to maximize 
the four-bond magnetization transfer, and no strong 
cross-peak between the H6 and HI' was observed in the 
corresponding 2D NOESY (D 2 0: 99.96%) spectrum, 
excluding the possibility of a syn x dihedral angle. 

Solution structures of UUU and 4"PT 

A total of 241 and 224 NOE restraints were used in the 
structure calculations of UUU and *4'\f, respectively 
(Table 1). Ten converged lowest-energy structures from 
each oligonucleotide were superimposed. An averaged 
pair-wise all-atom RMSD within the structure family of 
UUU is 1.15 A and that of I'** is 0.66 A (Figure 3a and 
b). In these structures, UUU and are both folded 

into stem-loop structures. To compare the structural 
effects of 4* modifications, the lowest-energy structures 
of UUU and were chosen to represent the 

NMR solution structures of the two oligonucleotides 
(Figure 3c and d). 

UUU and fff stem regions are folded into A- form 
duplexes containing a G'U wobble pair 

The stem regions (i.e. residues G1906-G1910 and C1920- 
C1924) of UUU and *W are folded into an A-form 
duplex conformation (Figure 3), similar to that of H69 
in X-ray crystal structures 1NKW (isolated large riboso- 
mal subunit from Deinococcus radiodurans with the same 
sequence and 4* modifications of H69 as E. coli) and 2I2T 
(large ribosomal subunit from E. coli assembled in 
complete ribosomes) (59,60). In the crystal structures of 
H69 from both bacteria and eukaryotes, G1907 and 
U1923 (E. coli numbering) form a wobble pair with two 
hydrogen bonds (G1907O6-U1923N3H and G1907N1H- 
U192302) (17,59-61), and this structural feature was also 
observed in the solution structures of UUU and 
(Supplementary Figure S5). The interproton distances 
between the imino proton of G1922 and the imino 
protons of G 1907 and U1923 (Supplementary Figure S5) 
support an A-form RNA base-stacking arrangement 
involving the G1907-U1923 wobble pair (44). 

*¥1911 and A 191 9 form a canonical Watson— Crick base 
pair with an additional intraresidue hydrogen bond 
involving V1911N1H 

Neighboring the stem-closing base pair (G1910-C1920), 
U1911 and A1919 in UUU do not form an optimal 
Watson-Crick base pair (Figure 4a). In the family of 10 
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Figure 2. 2D NOESY spectra of UUU (a) and *** (b) dissolved in H 2 0/D 2 0 (90/10%). Assignments of the base imino protons and cross-peaks of 
imino proton pairs are indicated. 



UUU NMR structures, the average distances between 
hydrogen-bond partner pairs UN3H-AN1 and U04- 
AN6H are both 2.5 A (2.5 ± 0.3 A and 2.5 ± 0.6 A, re- 
spectively), whereas the optimal hydrogen-bond distance 
is 1.8 A. These longer hydrogen-bond distances may reflect 
dynamics in the U191 1-A1919 pair or a paucity of NMR 
restraints. The sharp line widths for non-exchangeable res- 
onances at these positions suggest that the distorted U-A 
geometry is due to a lack of NMR restraints. However, it 
is worth noting that the H69 construct with U1911 is 
l.Okcal/mol less stable than a construct with *IT911 (37). 
In contrast, *1911 and A1919 in form a canonical 

Watson-Crick base pair (Figure 4b). The average dis- 
tances of vI>N3H-ANl and *02-AN6H are 1.8 ± 0.1 A 
and 1.7 ±0.2 A, respectively, and the corresponding 
three-atom geometry is much closer to a linear arrange- 
ment (170 ± 7° for *N3-MM3H-AN1 and 164 ± 6° for 
VI/02-AN6H-AN6). Formation of the *1911-A1919 base 
pair excludes solvent from the Watson-Crick edge of the 
two bases, leading to protection of these imino groups 
from CMCT and DMS chemical probing in the wild- 
type pseudouridylated 50S ribosomal subunit (40,62). 
Pseudouridine modification has also been shown to 
promote the formation of an optimal Watson-Crick 
base pair in tRNA Tyr from Bacillus subtilis to lock down 
the stem-loop transition and alter the conformation of the 
neighboring C*A + base pair on the loop side (55). 

The local structural rigidity of the stem-loop transition 
region of H69 may be further enhanced by potential 
water-mediated hydrogen-bonding interactions involving 
*1911N1H (36). Even though *1911N1H is exposed to 
the solvent (Supplementary Figure S6) and no partner is 



readily available within the direct hydrogen-bonding 
distance, the presence of a relatively sharp ^^UNIH 
resonance in the imino-proton spectrum suggests that 
this proton is protected from solvent exchange by a 
hydrogen bond, possibly mediated by a water molecule 
(2,36,58). In the family of 10 solution structures, 

the *1911N1H-02P distance is 4.1 ± 0.1 A and the Nl- 
N1H-02P angle is 115 ± 1° (Supplementary Figure S7a). 
Formation of a hydrogen bond involving ^1911N1H, 
together with base pairing of *IT911 and A1919, may 
explain the large chemical shift difference observed 
between the H6s of U1911 and ^1911 (Supplementary 
Figure S3). A similar local geometry has been observed 
in other pseudouridylated RNAs (6) and the crystal struc- 
tures of H69 in the isolated large subunit and 70S ribo- 
somes (PDB IDs 1NKW and 2I2T, Supplementary 
Figure S7b and c) (59,60). A thermodynamic stabilizing 
effect of * modification at 1911 has also been reported, in 
the context of a complete H69 hairpin (—1.0 kcal/mol) (37) 
and double-stranded constructs mimicking the stem 
region (residues from G1906 to A1912 and A1919- 
C1924 of H69, -1.1 kcal/mol) (38). These values agree 
with a thermodynamic contribution of —0.7 kcal/mol in 
conformational stability resulting from a water-mediated 
hydrogen bond to *N1H (36,58,63,64). 

W1915 and W1917 promote local base stacking 

Two of the three ^ modifications are clustered in the 3' 
half of the H69 loop, where distinct conformational dif- 
ferences are observed between the NMR structures of 
UUU and *W (Figures 4c and 5d). In residues 
4T915, A1916, vM917 and A1918 establish a continuous 
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Figure 3. NMR solution structure families (a and b) and the lowest 
energy structures (c and d) of UUU (a and c) and *** (b and d) are 
shown. In each of the structure families, 10 of the converged lowest 
energy structures were superimposed by 'alignment' with PyMol (The 
PyMol Molecular Graphics System, Schrodinger, LLC). The residues 
subjected to * modifications are labeled in (c and d). The backbone 
trace is shown as a gray ribbon, the base moieties are shown in ball and 
stick representation and the remaining are shown as a wire frame. 



base-stacking motif, whereas two disruptions in base 
stacking at the U1915-A1916 and U1917-A1918 steps 
are observed in the UUU RNA. These conformational 
differences are also consistent with the chemical-shift 
changes of the corresponding protons in this region 
when UUU and 4*4*4* are compared (Supplementary 
Figure S3) (54,65). In several other studies, 4*s were 
shown to have a propensity for base-stacking interactions 
(6,56,58). The step-wise distances of base mass centers of 
4*1915, A1916, 4*1917 and A1918 in 4*4*4* are 4~5A, 
within the distance range observed for base-stacking inter- 
actions in crystal structures (3.5~6A) (59,61). In the 
NMR structure of UUU, the distances between base 
mass o centers of U1915-A1916 (~7A) and U1917-A1918 
(~9A) steps are much larger. Because evidence for a 
QT-endo sugar pucker of U1915 and U1917 was 
revealed in the 2D DQF-COSY spectrum of the UUU 
construct (Supplementary Figure S4a), it is not surprising 
that these two residues can span much longer distances in 
UUU compared with the corresponding 4*s in the 
modified RNA (6,66). Given that the energy of base 



stacking is approximately proportional to 1/R 6 
(i.e. London dispersion interactions, in which R is the 
distance between the two base mass centers), the long dis- 
tances could decrease the stability of the base-stacking 
interactions between U 191 5-A 1916 and U 191 7-A 1918 in 
UUU (67). An additional 150 A 2 of solvent accessible area 
from the bases in the 4*1915 to A191 8 region of 4*4*4* is 
buried due to base stacking. There is no NMR evidence 
indicating involvement of the imino groups of 4*1915 and 
M>19 17 in hydrogen-bonding interactions. The conform- 
ational rigidity introduced by continuous base stacking 
promoted by 4*1915 and 4*1917 appears on first analysis 
to be inconsistent with previous thermodynamic studies 
(37,38), which show that 4*1915 and 4*1917 are slightly 
destabilizing. However, further analysis suggests that sta- 
bilization on the 3' part of the loop is compensated by 
destabilization on the 5' part of the loop. 

In previous studies, it was shown that pseudouri- 
dylations at 1915 and 1917 exerted a destabilizing effect 
of 0.8 kcal/mol on the global folding of H69 4*4*4* (38), 
with an individual thermodynamic contribution of 0.7 and 
0.3 kcal/mol from 4*1915 and 4*1917, respectively (37). 
The fact that 4* modifications are able to play opposing 
thermodynamic roles in the same RNA construct was un- 
anticipated. A rationale for these effects is that bending of 
the single-stranded region into a loop conformation is 
more demanding in the stacked environment. Goddard 
et al. reported that an extra 0.5 kcal/mol of enthalpy 
penalty is paid for each AA step in a DNA hairpin con- 
taining a polyA compared with polyT loop. This effect is 
attributed to base stacking of the purines. In contrast, in a 
DNA hairpin containing a polyT sequence, whose propen- 
sity to form base stacking is low, the thermodynamic 
penalty is purely entropic (68). The 4* modifications at 
1915 and 1917 appear to behave in a similar manner to 
the polyA hairpin loop, with Gibbs free energy changes 
(0.7 and 0.3 kcal/mol for 4*1915 and 4*1917) close to the 
reported value (0.5 kcal/mobbase stacking). This hypoth- 
esis is further supported by the fact that no imino protons 
of 4*1915 and 4*1917 participate in hydrogen-bonding 
interactions, as revealed in both the NMR spectra 
and the resulting structures, suggesting that the thermo- 
dynamic effects are dominated by base stacking. In 
contrast, in favor of hairpin folding of H69 4*4*4*, both 
N1H and N3H imino protons of 4*1911 are involved 
in hydrogen-bonding interactions that stabilize the 
4*191 1-A1919 base-pair platform, the terminal residues 
of the loop. 

A continuous base-stacking motif from 4*1915 to A19 18 
was observed in crystal structures of the eubacterial ribo- 
somal large subunit before (1NKW, D. radiodurans) and 
after (212T, E. coli) being assembled into complete ribo- 
somes (Figure 4e and f) (59,60). Because the sequence and 
4* modifications in H69 from D. radiodurans are the same 
as those in E. coli, the crystal structure of H69 from 
1NKW is used to represent the structure of H69 from 
E. coli before ribosome assembly, which is not currently 
available. Comparisons between each dinucleotide step in 
the H69 loop region of 1NKW and 2I2T indicate that, on 
assembly of complete ribosomes, residues in the 5' half 
undergo much larger conformational changes than 
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Figure 4. Local conformations of structural motifs involved in * modification(s) in NMR structures and crystal structures are shown. Unlike U1911 
and A1919 in UUU (a), *1911 and A1919 in *** (b) are positioned to form an optimal Watson-Crick base pair. In base moieties of residues 

from *1915 to A1918 form a continuous-stacking conformation (d), which is not observed in the same region of UUU (c). This continuous base 
stacking is also observed in the crystal structures of D. radiodurans (PDB ID 1NKW) (e) and E. coli (PDB ID 2I2T) (f) ribosomes. 



residues in the 3' half (black bars in Figure 5). Similar 
trends are observed when the NMR structure of 
and H69 crystal structure from 2I2T are compared 
(white bars in Figure 5). This observation is not 
surprising, as all three H69 RNAs (1NKW, 2I2T and 
contain the conserved *IT915 and <P1917 residues 
in the 3' half of the loop region, enhancing local base 
stacking. Comparison of the NMR structure of UUU 
and H69 crystal structure from 2I2T reveals large struc- 
tural differences in both the 5' and 3' halves of the loop 
region (gray bars in Figure 5), suggesting that more exten- 
sive conformational changes are required in H69 on 
ribosome assembly when *I* modifications are absent. 
A continuous base stacking in the 3' half of the H69 
loop promoted by ^ modifications is potentially able to 
provide consistent structural support, which may play an 
important role in maintenance of the global conform- 
ation during the highly dynamic translation processes 
(18,21-25). 

Effects of pseudouridine modifications on the 5' half of 
H69 loop region 

A general A-form RNA conformation is extended into the 
5' half (residues from A1912 to C1914) of the loop regions 
of UUU and (Figure 3). One major difference 

between the crystal structure 1NKW (isolated large 



ribosomal subunit from D. radiodurans) (59) and the 
NMR structures is that in 1NKW the base moiety of 
C1914 is stacked on *1915; in contrast in the NMR 
structures of UUU and >7.6A (the distance between 
A1913H8 and C1914H6) away from the base ring of 
A1913; in contrast in the NMR structures of UUU and 
the distances between the bases of A1913 and C1914 
are much closer (about 4~5A). The proximity of A1913 
and C1914 base moieties in space is confirmed by the 
medium or weak intensity cross-peaks between A1913H8 
and C1914H6 in the 2D NOESY spectra of UUU and 
However, no cross-peaks corresponding to 
C1914H6 and U/*1915H6 are observed in the same 
spectra, indicating a breakage of stacking between the 
bases of C1914 and U/M915 in the NMR structures 
(Figure 3), which is also observed in the crystal structure 
of a ribosome bound to the YaeJ protein (69). The differ- 
ent positioning of the C1914 base in NMR (UUU and 
*>I'*I') and crystal (1NKW and 2I2T) structures suggests 
that H69 assumes a slightly different conformation in 
solution, possibly reflecting the promotion of enhanced 
base stacking of C1914 and 4T915 under crystallization 
conditions. The difference in conditions of the NMR ex- 
periments and the crystallization studies (e.g. total salt 
concentration, Mg 2+ concentration, pH and solution 
versus crystal state) may also contribute to the differences 
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observed among the NMR structures and the crystal 
structure (Supplementary Figure S8). Chemical probing 
studies of whole ribosomes and large subunits also 
indicated that H69 changes conformation on varying the 
solution pH and Mg 2+ concentration (40,62). 




Dinucleotide step 

Figure 5. Conformational changes in each dinucleotide step on 
ribosome assembly are compared. Coordinates of the two neighboring 
residues in the NMR structures (UUU and and 1NKW were 

superimposed with corresponding residues in 2I2T and the RMSD for 
the two residues were calculated with all atoms. Dinucleotide steps were 
chosen to demonstrate the relative conformational changes within the 
two residues, when a propagating effect of coordinate translation 
resulted from other residues, e.g. stem duplex residues, is minimized. 



Residue C1914 is not the only nucleotide in the 5' half 
of H69 with conformational flexibility; several reports 
have indicated that A1913 flips out to form essential inter- 
actions with h44 of the 16S rRNA on assembly of the full 
ribosome to establish the B2a intersubunit bridge (16,59). 
Factors that can affect the structural dynamics of the 
5' half of the H69 loop region may play important biolo- 
gical roles. In summary, modifications at 1911, 1915 and 
1917 are capable of altering the conformational behavior 
of the 5' half of the H69 loop. Even though the global 
thermodynamic stability of is similar to that of 

UUU (AAG° 37 = -0.2kcal/mol) (70), this zero-sum 
result possibly elicits a conformational effect in the 
bridging region between M911 and the 3' half of the 
loop containing *1915 and ^1917. 



CONCLUSIONS 

Helix 69 is a 19-nt stem-loop structure in bacterial 23S 
rRNAs. The sequence of H69 is highly conserved within 
each domain of life, and its secondary structure is 
conserved throughout phylogeny (11). Post-transcriptional 
modification patterns in H69 also show conservation in all 
three domains of life (9,10). Even though * modifications 
in H69 are not essential for the survival and function of bac- 
teria (71), they confer significant advantages to the organism 
(35,72). Through comparisons of NMR solution structures 
of UUU and *I / <1'*I / , the modifications were shown to alter 
the structure, dynamics and thermodynamic behavior of 
the H69 loop (Figure 3). 





Figure 6. Structural comparisons of the loop regions of H69 UUU (left) and (right). *1915 and *1917 alter the loop folding by promoting 

local base-stacking interactions, and the conformational effects are locked within the loop region by a stabilized *1911-A1919 base pair, together 
with possible water-mediated hydrogen-bonding interactions involving *1911N1H. The conformational behavior of residues in the 5' half of H69 
loop region is affected by pseudouridylation. 
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Pseudouridine modifications in H69 are located in two 
structural regions. One pseudouridylation is at position 
1911, which is co-localized with A1919 at the H69 stem- 
loop junction. Formation of a canonical Watson-Crick 
base pair (*1911-A1919), as well as a possible water- 
mediated hydrogen-bonding interaction involving 
*1911N1H, helps maintain the stem duplex structure. 
These interactions allow conformational changes of H69, es- 
pecially those in the loop region, during the highly dynamic 
translation process (73). Conformational changes that occur 
in the H69 loop as a result of the presence of the * 

modifications at positions 1915 and 1917 are restricted to this 
region, while the stem structure is minimally perturbed 
(37,38,74). Two additional modifications, *1915 and 
*1917, occur in the 3' half of the loop region. Significant 
differences between the NMR structures of H69 UUU and 
are observed in this region. Without ^ modifications, 
only A1916 and U1917 participate in base stacking, and there 
are breaks in the interactions between the U1915-A1916 and 
U1917-A1918 steps. This local conformation of H69 has not 
been observed previously in crystal structures of the 
ribosome, as the RNAs in those structures contained the 
conserved *I* modifications. This 'unstacked' arrangement 
likely allows U1915, A1916 and U 191 7 to be more exposed 
to solvent and accessible to RluD, the enzyme responsible for 
converting U to * (34,75). After modifications are 
incorporated, the arrangement of residues <IT915, A1916, 
M>1917 and A1918 is 'stacked', in which a continuous base- 
stacking motif is established. Re-organization of the local 
conformation on pseudouridylation not only reinforces the 
structural rigidity in the 3' half of H69 WI* loop (76), but 
also evokes a structural effect that ripples upstream into the 5' 
half of the loop region to alter local conformational 
dynamics, which plays an important role in ribosome 
assembly and maintenance of the ribosomal structural integ- 
rity during the translation process. 

In summary, exposure of loop residues U1915, A1916 
and U1917 is likely important for substrate recognition 
and catalysis by RluD. After incorporation of ^ modifi- 
cations at positions 1911, 1915 and 1917, these nucleotides 
appear to work synergistically in modulating the conform- 
ation and structural dynamics of the loop (Figure 6), 
which plays a role in formation and maintenance of 
bridge B2a during various stages of translation. It has 
been reported in studies of telomerase (6), spliceosomal 
(7) and transfer (8) RNAs that * helps pre-structure 
local RNA motifs for subsequent participation in down- 
stream events; therefore, this work on H69 helps complete 
the spectrum of roles for * in all major biologically 
relevant RNA species. The unique structural properties 
of * and its ability to fine-tune the structure of function- 
ally important RNAs make these sites attractive targets 
for drug design, and the structures obtained in this work 
can be applied in the development of novel inhibitors of 
bacterial ribosome function. 
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